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In a receiver of a digital radio communication system a combined channel estimate Oup) is fonned by averaging a long channel 
estimate (hk4 with a short channel estimate (hk^). 



BNSOOCID: <WO_9613910A1J^ 







FOR THE PURPOSES OF INFORMATION ONLY 






Codes used to identify States paity to the PCT on the from pages of pamphlets publishing international 


applications under the PCT. 










AT 


Austria 


GB 


United Kingdom 


MR 


Mauritania 


AU 


AutcnlU 


GE 


Geofgia 


MW 


Malawi 


BB 


Bartwdos 


GN 


Guinea 


NE 


Niger 


BE 


Belgium 


GR 


Greece 


NL 


Netherlands 


BF 


Bufkina Faso 


HU 


Hungaiy 


NO 


Norway 


BG 


Bulgaria 


IE 


Ireland 


NZ 


New Zealand 


BJ 


Benin 


IT 


Italy 


PL 


Poland 


BR 


Brazil 


JP 


Japan 


PT 


Portugal 


BY 


Belarus 


K£ 


Kenya 


RO 


Romania 


CA 


Canada 


KG 


Kyrgysxan 


RU 


Russian Federation 


CF 


Centra] African Republic 


KP 


Democratic Pfeopk'i Republic 


SD 


Sudan 


CG 


Congo 




of Korea 


SE 


Sweden 


CH 


Swiucrland 


KR 


Republic of Korea 


SI 


Slovenia 


CI 


COie d* I voire 


KZ 


Kazakhstan 


SK 


Slovikia 


CM 


Cameroon 


U 


Liechtenstein 


SN 


Senegal 


CN 


China 


LK 


Sri Lanka 


TD 


Chad 


CS 


Czechoslovakia 


LU 


Luxembourg 


TC 


Togo 


CZ 


Czech Republic 


LV 


Latvia 


TJ 


Tajikistan 


DE 


Germany 


MC 


Monaco 


XT 


Trinklad and Tobago 


DK 


Denmark 


MD 


Republic of Moldova 


UA 


Ukraine 


ES 


Spain 


MG 


Madagascar 


US 


United States of America 


FI 


Fmland 


ML 


Mali 


UZ 


Uzbekistan 


FR 


Frmoe 


MN 


Mongolia 


VN 


Viet Nam 


GA 


GabOD 











BNSOOCID: <WO 9613910A1J_> 



wo 96/13910 PCT/SE95/01275 

1 

METHOD AND APPARATUS FOR CHANNEL ESTIMATION 
TECHNICAL FIELD 

The present invention relates to a method and apparatus for 
forming a channel estimate in a digital radio communication 
5 system. 

BACKGROUND OF THE INVENTION 

In TDMA radio communication systems (TDMA = Time Division 
Multiple Access) and other frame based communication systems 
information is transmitted on a channel in the form of signal 

10 frames or bursts . In order to synchronize the receiver to these 

signal frames each signal frame comprises a known synchronization 
word in a predetermined position within each signal frame. In for 
instance the European GSM system for mobile telephony this 
synchronization word is 26 bits long. When the receiver expects 

15 a new signal frame from the transmitter, a training sequence that 

is identical to the transmitted synchronization word is generated 
by a training sequence generator in the receiver. The received 
signals are compared to the locally generated training sequence, 
and when the best possible correlation is obtained between this 

20 sequence and the received signals, synchronization is considered 

to exist between the locally generated and the received signal . 

In addition to the synchronization itself the training sequence 
is also used for channel estimation. Since the radio channel 
often is subjected to multiple path propagation the receiver 

25 comprises some sort of equalizer to eliminate this .phenomenon . 

The equalizer requires a time limited estimate of the impulse 
response of the channel . This impulse response can be obtained 
from the correlation signal. Forney [l] and Ungerboeck [2] 
describe two different algorithms that, given the channel impulse 

30 response and Gaussian channel noise with }cnown correlation 

proDerties, determine the most likely sent sequence. Both 
algorithms will work properly, but with degraded performance, if 
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an approximate estimate of the channel impulse response is used 
instead of the true impulse response, or if the noise is non- 
Gaussian. The equalizer makes use of the channel estimate to 
initialize and update e.g. filter taps. An example is the Maximum 
5 Likelihood Sequence Estimation (MLSE) detector, where an FIR- 

filter is used directly as a model of the channel [2]. Another 
example is decision feedback equalizers (DFE) [4], where the 
filter taps in both forward and feed-back are calculated from the 
channel estimate. 

10 When the channel is estimated from a received synchronization 

word the estimate will contain noise, since the received data is 
noisy and the training sequence is of finite length. Even a 
channel estimate that is updated continuously will be noisy. The 
noise content in the channel estimate will be very high in e.g. 

15 fading dips, since the signal level in that case is low compared 

to the noise or interference level, giving a low signal-to-noise 
ratio in the received data. The synchronization process is also 
disturbed by receiver noise. An incorrect and unstable synch- 
ronization gives an incorrect channel estimate as a secondary 

20 effect, thereby causing a substantial loss in receiver performan- 

ce. The noisy channel estimate gives incorrect settings for the 
filter taps in the detector or equalizer, e.g. the MLSE will have 
an incorrect FIR channel model. This causes degraded demodulation 
performance in a way similar to the degradation from noisy 

25 received samples. The effect is that the interference or noise in 

the received samples has a double impact - first the channel 
model is impaired, then the incorrectly set equalizer will have 
to equalize and demodulate noisy samples. 

A common feature of the described prior art is that the channel 
30 estimate is used directly in the detection or equalization 

procedure without considering the noise content of the estimate. 
However, for example the MLSE is optimized in the sense that it 
gives the most likely demodulated symbol sequence only if the 
receiver noise is additive Gaussian and the channel estimate is 
35 exact. It is not optimal if there is noise in the channel 
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3 



10 



estitnate . 

A .«chod Of reducing the influence of noisy taps in the channel 
es'^te I proposed in US patent S.25X.233 ..abedz et al, . There 
T^ TslsJJro delete taps in the channel estimate that are 
: lc„ a Lrtain threshold value, thereby «f 
contribution from noisy taps with low content of useful sxgnal . 
Totally eliodnating some channel estimate taps may. however, 
™ Vital information, since it is very difficult to dis- 
tinguish between useful signal energy and noise energy xn a tap. 

EP-A-O 535 403 describes a method in which a channel estimate of 
a current burst is combined with a channel estimate of a previous 
burst. Both estimates have the same number of taps. 

SUMMARV OF THE INVENTION 

Thus, an Object of the present invention is to provide a method 
and apparatus for forming an in^roved channel estimate :Ln a 
digital radio communication system. 

The method of the invention is characterized by the features of 

claim 1 . 

The apparatus of the invention is characterized by the features 

20 of claim 9. 

The present invention is based on the observation that a long 
channel estimate (an estimate having many taps, will contain more 
noise than a short estimate .having fewer taps, . This xs shown xn 
Che APPENDIX. Where the tap noise is given by equatxon (9,^ In 
this equation N represents the total number of samples used xn 
the calculation of the channel estimate and M represents the 
number of taps in the channel estimate. Increasing the number of 
taps M in the channel estimate will therefore also xncrease the 
tao noise Furthermore, the total noise contrxbutxon xn the 
deLdulatxon process in the detector or e<pializer is also 



15 
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proportional to the number of taps in the channel estimate. From 
this point of view it is therefore desirable to have a channel 
estimate that is as short as possible. 

On the other hand, a short channel estimate may, in a case where 
5 there is much time dispersion, cause the receiver to exclude 

received energy outside of the span of the estimate, thereby 
disregarding important information present in the signal. The 
short channel estimate does, however, include the part of the 
signal with most of the energy content. 

10 The solution to this problem is, according to the present 

invention, to combine one long channel estimate with at least one 
shorter channel estimate- This combines the properties of both: 

1. It contains reliable (not so noisy) information on the part 
of the received signal that has most of the energy. 

15 2 . It also contains information on the signal spread caused by 

time dispersion over a larger time interval . 

BRIEF DESCRIPTION OF THE DRAWINGS 

The invention, together with further objects and advantages 
thereof, may best be understood by making reference to the 
20 following description taken together with the accompanying 

drawings, in which: 

FIGURE 1 is a block diagram of an apparatus in accordance 

with the present invention; 

FIGURE 2 is a flow chart illustrating the method in accor- 

25 dance with the present invention; 

FIGURE 3 illustrates how channel estimates are combined in 

one embodiment of the present invention; 
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FIGURE 4 illustrates how channel estimates are combined in 

another embodiment in accordance with the present 
invention; 

FIGURE 5 illustrates how channel estimates are combined in 

a further embodiment in accordance with the present 
invention; and 

FIGURE 6 illustrates a simplified calculation of the ampli- 

tude of a complex number, which may be used in a 
preferred embodiment of the invention. 

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS 

The apparatus and method of the present invention will now be 
described with reference to Figs, l and 2. 

In Fig. 1 an A/D-converter receives analog samples b(nT.) and 
converts these samples to a sequence of digital samples b„. These 

15 digital samples b„ are forwarded to a Maximum Likelihood Sequence 

Estimator 12, which outputs a sequence of detected symbols u„. 
The digital sequence b„ is also forwarded to a correlator 16 that 
correlates sequence b„ with a locally generated training sequence 
received from training sequence generator 17 . The correlation 

20 values c, from correlator 16 are used to synchronize the burst 

(step 110 in Fig. 2). This synchronization step will be further 
described below. 

Methods for burst synchronization will be described below with 
reference to the European GSM- system. In this system a synch- 

25 ronization word comprises 26 bits. The 16 central bits in this 

word have good correlation properties when correlated with the 
entire synchronization word, that is a maximum correlation = 16 
in the central position and a correlation of 0 in the remaining 
ten positions (C{k) = [0 0 0 0 0 16 0 0 0 0 01 ) . These 16 central 

30 bits are generated as a training sequence in a training sequence 

generator in the receiver. This training sequence is used for 
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forming for instance 11 correlation values Ci with the received 
signal frame. According to one burst synchronization method the 
final synchronization position is chosen by comparing mutually 
displaced windows, each containing 5 correlation values, with 
5 respect to energy contents, and choosing the time position of the 

window with maximum energy as the synchronization position. 
Another burst synchronization method is described in EP-A-0 551 
803. 

Since two channel estimates will be combined, burst synch- 
10 ronization is performed for both the short and the long estimate. 

The synchronization step is performed in synchronizers 18 and 19, 
respectively. (Since the estimates are of different length they 
will not necessarily synchronize to the same burst position.) 

The synchronization positions are forwarded to channel estimators 
15 20, 2 2 for the long and short channel estimate, respectively. 

These estimators calculate channel estimates around the respec- 
tive synchronization positions, as will be further described 
below . 

In the APPENDIX it has been shown that the tap noise may be 
20 estimated in accordance with formula (9) . This formula indicates 

that the tap noise will be reduced by using as many samples b„ as 
possible (increasing N) . For this reason the long estimate h^'' is 
recalculated by using as many of the 26 samples of the training 
sequence u^ as possible and the received synchronization word b„. 
25 Thus, channel estimator 20 will calculate five correlation values 

(in GSM) from N-M+l= 26-5+1=22 samples (all 26 samples are 
actually used since 5 correlation values each based on 22 samples 
and shifted one sample are formed) . These calculations are 
performed in the same way as in correlator 16, but since we now 
30 know the synchronization position the entire training sequence 

may be used to form the five taps of the long channel estimate 
hk**. The described procedure corresponds to step 120 in Fig. 2. 

The short channel estimate h/ is formed in a similar way in 
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Channel estimator 22. However, in this case fewer than five taps 
have to be calculated (M<5) . In a preferred embodiment of the 
present invention the short channel estimate only comprises one 
tap, which means that N-M+l«26-l*l=26 samples may be used for 
5 chamiel estimation. Thus, in this case the entire training 

sequence is used to form a single correlation value with a 
significantly reduced tap noise. The described procedure 
corresponds to step 130 in Fig. 2. 

The calculated channel estimates h^*" and h^» are combined in an 
10 averaging circuit 24 (corresponds to step 150 in Fig. 2). An 

example of this averaging process is illustrated in Fig. 3 (to 
simplify the illustration only the amplitudes of the estimates 
are shown) . In this embodiment the estimates are combined with 
their original burst synchronization (since the illustrated short 
15 estimate only has one tap it is assumed to have taps with zero 

values in all other time positions) . This may result in a 
situation where the maximum values of the two estimates are not 
in the same time position, as indicated in Fig. 3. 

Since the short and long channel estimates have been separately 
20 synchronized with the burst, it may be preferable to first 

mutually synchronize the calculated channel estimates. This is 
indicated by line 26 in Fig. l. where channel estimator 20 
informs synchronizer 18 of the proper synchronization position to 
be used for the short channel estimate. This corresponds to step 
25 140 in fig- 2. 

There are different ways to obtain mutual synchronization between 
the short and long channel estimates. One method is illustrated 
by Fig. 4- According to this method the maximum tap of the short 
estimate h^' is synchronized with the position of the maximum tap 
30 of the long estimate h^*-. This embodiment implies that the taps 

of the short estimate may be shifted, as indicated by the dashed 
tap in Fig. 4. This dashed tap represents the burst synchroniza- 
tion position of the short estimate. The tap will be recalculated 
in estimator 22 and then shifted to the position of the solid 
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line tap estimate h^^ in Fig. 4. This estimate synchronization is 
reasonable, .since the maximum taps of the long and short 
estimates usually will have coinciding positions. 

Another mutual channel estimate synchronization method is 
5 illustrated in Fig. 5. Here the burst synchronization position of 

the short estimate (containing only one tap in the example) h^* 
is indicated by the dashed tap. However, this tap will not be 
recalculated. Instead the tap that corresponds to the position of 
the maximum tap of the long estimate h^*- will be recalculated and 
10 averaged with the maximum tap of the long estimate. As in the 

previous embodiment this estimate synchronization is reasonable, 
since the maximum taps of the long and short estimates usually 
will have coinciding positions . 

After mutual synchronization of the short channel estimate h^* and 
15 the long channel estimate h^*- these estimates are combined in an 

averaging circuit 24. In a simple, presently preferred em- 
bodiment, this averaging circuit 24 performs sitt5>le averaging of 
corresponding taps of the long and short estimates, as illustra- 
ted in Figs . 3 , 4 and 5 . In a more elaborate embodiment a 
20 weighted average may be formed by calculating weighting factors 

from reliability measures (noise measures) for each of the two 
channel estimates. The combined channel estimate hk*= is forwarded 
to maximum lilcelihood sequence estimator 12 for adjusting the 
filter coefficients thereof. 

25 In case the short chaimel estimate contains only one tap, the 

above described process may be simplified by noting that the 
maximum correlation value c, may be used directly to represent 
this tap. Thus, if a reduction of calculation complexity is 
desired, it may be sufficient to provide the maximum value of Cj 

30 as the short estimate (no recalculation of the short estimate) . 

In such an embodiment the burst synchronization position of this 
estimate may be used as a reference position for the largest tap 
of the long estimate. 
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elaborate entooditnent the energy of the maximum tap of 
in a more el^ja^e e ^^^^^ ^^^^^ ^^^^ 

rtil't! Xf^rmLilm'tap is very .ominatine this indicates 
r^IHhere xs very little time dispersion. In this case a sxngle 
odtrof the^hannel is a very good model and the single tap 
5 tap model of the cha - ^^^i^e training 

short estimate xs recalcu .dispersion is large the 

seauence- On the other hand, xf tne txme ais.F 

sequence spread out and the maximum tap wxll not 

channel estxmate wxll ne spiea^a « 

channex , ' v,^<= .-ase a single tap model is not so good and 

e^Ire .raining science. Thus, in this case the sxngle tap fro™ 
the hurst synchronization is used as the short estimate. 

Mother Simplification that may be performed is to consider only 
Te amplitude of the short estimate and disregard its phase <.n 
case the short estimate contains only one tap) . This non-coherent 
co^ination does not in^ly a significant loss in performance 
sl^e the phase of the two largest taps of both estimates usually 
are close to each other. 

A further simplification may be performed by calculating the 
.0 U Of the short estimate by adding the ^^^^^ - 

its real and imaginary parts, respectively, as xllustrated in 
«g 7 in this approximation z is approximated by |x Myt , where 
X and y are the real and imaginary part, respectively, 
in tL illustrated embodiments two channel estimates have been 
co,^ined. However, it is obvious that it is possible to corr^ne 
25 comDiii« feasible embodiment would be a 

«,«T-o ^>^an two estxmates . Tnus, a j.e«=»j-"-i- 

7Z::Z1 Of three channel estimates of ^ 
is also appreciated that the present xnventxon xs applxcable 
other systems than the described GSM system. 

^ ™»^hods described above are simple ways to improve receiver 
Te flr^" Tthout dramatic increase in comp-xity^ The 
Improvement in receiver performance xs approxxmately 1 dB 
some propagation conditions in the GSM case. 
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It will be understood by those skilled in the art that various 
modifications and changes may be made to the present invention 
without departure from the spirit and scope thereof, which- is 
defined by the appended claims- 
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APPENDIX 



^rn.arec. e«;r i m^t^QT^ of t-»e channel 

The model used for t.he channel is a FIR filter with filter taps 
{h,} modelling the propagation and an additive white Gaussxan 
noise (AWGN) source e, modelling the noise. It is sxmxlar to 
Haykin's linear regression model (Simon Haykin. "Adaptive Fxlter 
Theory". Prentice -Hall . 1986, pp 307-316) . The expression for the 
received discrete signal samples is: 



10 



15 



20 



where u, is the input signal to the channel and M is the length 
of the channel estimate (filter taps {hj). Thus, the radio 
Channel impulse response is limited to M samples of h,. In GSM 
for example M=5 . The impulse response is an unknown parameter 
estimated from the training sequence. In the GSM example the 
training sequence is N«26 symbols long. A least squares method 
(see Haykin [3]) is used to estimate the M tap model (5 taps xn 
GSM) The input data from the training sequence is windowed usxng 
the covariance method. Thus, there are N-M.l (22 in GSM) received 
samples to use for estimation. The input is arranged as a matrxx 
A and the received samples b. error e and impulse response h are 
arranged as vectors as follows: 



•M-1 



U 



"iV-X 



... U„. 



(2) 
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b = [jbo 2>x ... b„.„f (3) 
d = [eo ... ej,.„f <4) 

i = [A, A, ... Vi]"" 

where T represents transpose. In vector form the channel model 
may be expressed (for the training sequence) as: 

jb = Ah + e <6) 
The least squares estimate £ of h is (see Haykin) 



where *=A"A is the deterministic correlation matrix of the 
training sequence. Here H represents Hermitean transpose. 

The variance of each element of the channel estimate hi is called 
Ni, and depends on the noise variance and on the training 
sequence used. Ni is the ith diagonal element of a^*"S where ^ is 
the deterministic correlation matrix defined above. 



10 If the training sequence is close to white (which is the case in 

e.g. GSM), the correlation matrix will be: 

N-M-^ 1 



where 1 is the identity matrix. The conclusion is that the tap 
noise is: 

N. = 2 (9) 
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CLAIMS 



15 



1. A method of forming a channel estimate in a digital radio 
communication system, characterized by: . ^ ^ 

forming a first channel estimate having a predetermined number 

of taps; 

forming at least one further channel estimate havxng fewer 
taps than said first channel estimate; and 

combining said first and said further channel estimates for 
forming a combined channel estimate. 

2 The method of claim 1, characterized by averaging correspon- 
ding taps of said first and said further channel estimates, saxd 
further channel estimate (s) being padded with zero values in tap 
positions that correspond to taps of said first channel estimate 
that have no correspondence in said further channel estimate (s) . 

3 The method of claim 2 , characterized by considering only tap 
amplitude and disregarding tap phase of taps of said further 
channel estimate (s) in said averaging. 

4 The method of claim 3, characterized by approximating the 
amplitude of each tap of said further channel estimate (s) wxth 
the sum of the amplitudes of its real and imaginary parts, 

respectively - 

5. The method of claim 2, 3 or 4, characterized by a single 
further channel estimate having a single tap. 

6 The method of claim 5, characterized by said single tap of 
said further channel estimate being averaged with a tap of 
maximum amplitude in said first channel estimate. 

7 The method of claim 5, characterized by averaging a tap of 
maximum an«)litude of said first channel estimate with the single 
tap of said further channel estimate, as estimated xn the same 

30 time position. 



20 



25 
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8. The method of claim 2. 3, 4, 5, € or 7, characterized by 
weighted averaging of said first and said further channel 
estimates - 

5 9. An apparatus for forming a channel estimate in a digital radio 

communication system, characterized by: 

first means (20) for forming a first channel estimate having 
a first predetermined number of taps; 

second means (22) for forming at least one further channel 
10 estimate having fewer taps than said first channel estimate; and 

means (24) for combining said first and said further channel 
estimate (s) for forming a combined channel estimate, 

10. The apparatus of claim 9, characterized by averaging means 
(24) for averaging corresponding taps of said first and said 

15 further chamnel estimate(s), said further channel estimate{s) 

being padded with zero values in tap positions that correspond to 
taps of said first channel estimate that have no correspondence 
in said further channel estimate (s) . 

11. The apparatus of claim 10, characterized by said averaging 
20 mesois (24) considering only tap amplitude and disregarding tap 

phase of taps of said further chamnel estimate (s) in said 
averaging . 

12. The apparatus of claim 11, characterized by said averaging 
means (24) approximating the amplitude of each tap of said 

25 further chcunjiel estimate (s) with the sum of the amplitudes of its 

real and imaginary parts, respectively. 

13. The apparatus of claim 10, 11 or 12, characterized by a 
single further channel estimate having a single tap. 

14. The apparatus of claim 13, characterized by said first means 
30 (20) synchronizing said single tap of said further channel 

estimate with a tap of msiximum amplitude in said first channel 
estimate . 
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15 The apparatus of claim 13 . characterized by said first means 
(20) synchronizing a tap of maximum amplitude of said first 
channel estimate with the single tap of said further channel 
estimate, as estimated in the same time position. 
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